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The Riemann hypothesis states that all nontrivial zeros of the zeta function lie on the critical 
line Re(s) = 1/2. Hilbert and Polya suggested a possible approach to prove it, based on spectral 
theory. Within this context, some authors formulated the question: is there a quantum mechanical 
system related to the sequence of prime numbers? In this Letter we assume that there is a class 
of hypothetical physical systems described by self-adjoint operators with countable infinite number 
of degrees of freedom with spectra given by the sequence of primes numbers. We prove a no-go 
theorem. We show that the generating functional of connected Schwinger functions of such theories 
cannot be constructed. 
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Prime numbers occur in a very irregular way within 
the sequences of natural numbers. In particular, the dis- 
tribution of prime numbers exhibits a local irregularity 
but a global regularity. The best result that we have 
concerning their global distribution is the prime number 
theorem: the number of primes 7r(x) not exceeding x is 
asymptotic to xj In x [lHy]. 

The Riemann hypothesis claims that all nontrivial ze- 
ros of the Riemann zeta function Q(s) lie on the critical 
line Re(s) = 1/2 Q- Since in the region of the com- 
plex plane where it converges absolutely and uniformly 
there is a representation for £(s) in terms of the prod- 
uct of all prime numbers, there is a connection between 
these nontrivial zeros of the f (s) function and the dis- 
tribution of prime numbers in the sequence of integers. 
One way to connect the distribution of primes and the 
non-trivial zeros of the Riemann zeta function is given 
by the integral logarithm function Li(a;). The assertion 
that ir(x) = Li(x) + 0(y/x\nx) for x — > oo is equivalent 
to the Riemann hypothesis 0. 

Hilbert and Polya suggested that there might be a 
spectral interpretation for the non-trivial zeros of the 
Riemann zeta function. The nontrivial zeros could be 
the eigenvalues of a self-adjoint linear operator in an ap- 
propriate Hilbert space. The existence of such operator 
leads us to the proof of the Riemann hypothesis. The 
spectral nature of the nontrivial zeros emerges studying 
the local statistical fluctuations of high zeros. These fluc- 
tuations are universal and follow the law of fluctuation 
of the eigenvalue distribution for matrices in the Gaus- 
sian unitary ensemble a result that prompted many 
authors to consider the Riemann hypothesis in the light 
of random matrix theory. 

Since the main point of Riemann's original paper is 
that the two sequences, of prime numbers on one hand, 
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and of Riemann zeros on the other hand, are in duality, 
any statement about one of these two sequences must 
have a translation in terms of the other. Therefore, based 
on the Hilbert-Polya conjecture some authors formulated 
the following question: is there a quantum mechanical 
system related to the prime numbers [lOl-QjJ? The idea 
is to assume that there is a trace formula for the prime 
numbers and therefore a Hamiltonian with a prime num- 
bers spectrum. For a nice discussion concerning physics 
of the Riemann hypothesis, see flU ]. 

Quantum field theory is the formalism where the prob- 
abilistic interpretation of quantum mechanics and the 
special theory of relativity were gathered to account for 
a plethora of phenomena that cannot be described by 
quantum mechanics of systems with finitely many de- 
grees of freedom. Euclidean methods lead to mathe- 
matically well-defined functional integral representation 
of the solutions of field theory models. Spectral prop- 
erties of a Hamiltonian and momentum operators of a 
field theory can be deduced from cluster properties of 
moment of the measure that determines the functional 
integral. Furthermore, boson field theory in one di- 
mension (time) reduces to quantum mechanics formu- 
lated in the Heisenberg picture. That is, any quantum 
mechanical system with n degrees of freedom derives 
from some n-componcnt scalar field model in one dimen- 
sion. A stochastic process canonically associated with a 
Hamiltonian Hq + V corresponds to Euclidean field self- 
interaction via the potential V 14 1 . 

Based on above discussions, one might be tempted to 
inquire whether quantum field theory techniques can be 
used to exploit questions related to number theory In 
such a direction, connections between field theory and 
prime numbers have already been investigated by many 
authors The aim of this Letter is to investigate 

the possibility of existence of a quantum field possessing 
a spectrum given by the sequence of prime numbers. For 
simplicity, we consider a free scalar field in Euclidean 
space, but the results derived here can be extended to 
higher spin fields. We assume that there is a hypothct- 
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ical class of physical systems described by self-adjoint 
operators associated with bosonic fields with a spectrum 
given by the sequence of primes numbers. We show that 
the functional integrals associated with such theories are 
ill-defined objects. Here we use h = c = 1. 

The Riemann zeta function £(s) is a complex variable 
function of s, i.e., s = a + it with a, r € 1Z, which is 
defined as the analytic continuation of a Dirichlct series 
to the whole complex plane. Its unique singularity is 
the point s = 1 at which it has a simple pole with unit 
residue. In the half-plane Re(s) > 1 it is defined by the 
absolutely convergent series 
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(i) 



In the complex plane for Re(s) > 1 the product of all 
prime numbers appears as a representation of the Rie- 
mann zeta function: 



(2) 



The product given by Eq. ((2]) is called the Euler prod- 
uct. This is an analytic form of the fundamental theorem 
of arithmetic. Since a convergent infinite product never 
vanishes, from Eq. ([2]) we obtain that C(s) has no zeros 
for Re(s) > 1. The proof that the zero-free region in- 
cludes Re(s) = 1, i.e., £(1 + it) ^ is another way to 
formulate the prime number theorem. From the func- 
tional equation 



K -^r(i^i)c(i- s ), (3) 



valid for s € C\{0, 1} it is possible to show that the region 
Re(s) < is also zero-free with respect to the non-trivial 
zeros. It is often more convenient to write it as 



C(s) 



r(f + 1) (s l) 



(4) 



The function £(s) is analytic everywhere on the complex 
plane and the poles of the Gamma function reproduce all 
the zeros of £(s) on the real axis at a = — 2k, k = 1, 2, .... 
Thus £(s) contains only the complex zeroes of ((s). There 
is no available proof that there are no other nontrivial 
complex zeros away from a = 1/2. 

Now let us construct a field theory satisfying the con- 
dition of a prime numbers spectrum. Consider a free 
neutral scalar field defined in a c?-dimcnsional Minkowski 
spacetime. The Euclidean field theory can be obtained 
by analytic continuation of the n-point Wightman func- 
tions to imaginary time. Let us assume a compact Eu- 
clidean space with or without boundary. If the compact 
Euclidean space has a smooth boundary, one imposes 
suitable boundary conditions on the fields. Let us sup- 
pose that there exists an elliptic, self-adjoint differential 
operator D acting on scalar functions on the Euclidean 



space. The usual example is D = (— A + m\ ), where A 
is the <i-dimensional Laplacian. The kernel K(niQ; x — y) 
is defined by 



K(m Q ; x-y) = (-A + ml ) 6 d (x - y). 



(5) 



The Euclidean generating functional Z[h] is formally de- 
fined by the following functional integral: 



Z[h] = J [dip] exp ^-S + J d d a 



h(x)ip(3 



(6) 



where the action that usually describes a free scalar field 
is 

S (p) = / d d xd d yp(x)K(m ;x-y)ip(y). (7) 



In Eq. (|6|), [dep] is a translational invariant measure, for- 
mally given by [dip] = Y\ x dip(x). The term m§ is the 
(bare) mass squared of the model. Finally, h(x) is a 
smooth function introduced to generate the Schwingcr 
functions of the theory. One can define the functional 
W[h] = \nZ[h] which generates the connected Schwingcr 
functions. In order to obtain a well-defined object, we 
need to regularize a determinant associated with the op- 
erator D, since W[0] = — l/21ndet D. A similar situation 
arises for the case of self-interacting scalar fields when one 
calculates the one-loop effective action. For both cases, 
a finite result can be achieved by defining the spectral 
zeta function associated with the corresponding elliptic 
operator (l9M2l|. 

The standard technique of the spectral theory of el- 
liptic operators implies the existence of a complete or- 
thonormal set {/fe} fc=1 such that the eigenvalues obey: 



< Ai < A 2 < 



< At 



when k — > oo. In 



the basis {fk} the operator D is represented by an in- 
finite diagonal matrix D = diag (Ai, A2, ...). Therefore 
the generic operator D satisfies Df n {x) = \ n f n {x). The 
spectral zeta function associated with the operator D is 
defined as 
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Cd(s)=X!t7' Re(s) > s , 



for some sq. Formally we have 



— C,d{s)\ s =q 
as 



lndetD. 



(8) 



(9) 



From equation ([9]), one notes that the spectral zeta func- 
tion has to be analytically continued beyond the domain 
of convergence of the series. In other words, it is neces- 
sary to perform an analytic continuation of Cd(s) from 
an open connected set of points for sufficient large pos- 
itive Re(s) into the whole complex plane. In particular, 
the spectral zeta function must be defined at s = 0. This 
method can also be employed in the case of interacting 
fields in curved space, at least in a ultrastatic space-time. 
In this case the one-loop effective action is computed by 
evaluating a determinant, as mentioned above. Since the 
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argument of the logarithm must be strictly dimension- 
less, one may have scaling properties: 

5^M- 2 o( s )l*=o = ^ hi fJ?( D \s=Q + ^Cd|s=o, (10) 

where /i is an arbitrary parameter with dimensions of 
mass. 

Let us assume that there is a hypothetical operator 
D acting on scalar functions in the Euclidean manifold 
with prime numbers as its spectrum. We are using the 
idea of relating a numerical sequence with the spectrum 
of a differential operator. In this case, the spectral zeta 
function associated with the hypothetical operator D is 
the so-called prime zeta function P(s), s = a + it, for 
real a and t, defined as 



p(a) 



Re(s) > 1. 



(11) 



The above summation is performed over the sequence of 
all primes [H, [23[. Since ln(l — 1/p) is of order — 1/p 
and Yip (1 — P^ 1 ) = 0, the sum of the reciprocals of the 
prime numbers diverges, but the series converges abso- 
lutely when a > 1. From the above discussions, in or- 
der to regularize the generating functional of connected 
Schwinger functions, we must analytically extend P(s) to 
s = 0. Using the generalization for the Euler result given 
by Eq. ([2]) and the definition of the prime zeta function 
we have 



oo 

lnC(s) =J2~ P ( rs )' Re ( s ) > L 



Employing the Mobius inversion formula, one has 



(12) 



P(*) = E^ lnC(fcs) ' Rc ( s )> 1 ' ( 13 ) 

fc=l 

where the Mobius function /i(n) is defined as [24[ 

1 if n is a square-free positive integer 
with an even number of prime factors, 
fj,(n) = ^ —1 if n is a square- free positive integer 

with an odd number of prime factors, 
if n is not square-free, 

where /i(l) = 1. From Eq. ([13]) one notices that the 
analytic continuation of the prime zeta function can be 
obtained from the analytic continuation of the Riemann 
zeta function. Introducing the Jacobi theta function 



J2 e-" 2 ™, x>0 



(14) 



n— — oo 



and using the usual definition of the Gamma function, 
we have 



'I 2 ) TT-iCOO 



s(s 



(15) 



where i?i(x) = \(&{x) — 1). The integral that appears in 
the above equation is convergent for all values of s and 
therefore this equation gives the analytic continuation 
of the Dirichlet series defined by Eq. (JTJ) to the whole 
complex s-plane. The only singularity is the pole at s = 
1, since the pole at s — is canceled by the pole of the 
Gamma function T(|). Therefore the analyic extension 
is valid for s G C \ {1}. 

The analytic structure of P(s) for Re(s) < 1 is given by 
the poles and nontrivial zeros of the Riemann zeta func- 
tion. First, due to the polar structure of the Riemann 
zeta function, in Eq. (| 13[) we get that s = ^ is a singular 
point for all square free positive integers k. This sequence 
limits to s = 0. All points on the line Re(s) = are limit 
points of the poles of P(s), so that the line Re(s) = is 
a natural boundary of P(s). As discussed in [22(, there 
is a clustering of singular points along the imaginary axis 
from the non-trivial zeros of the Riemann zeta function. 
Hence the prime zeta function has an analytic continua- 
tion only in the strip < a < 1. As discussed, the only 
possible singularities of P(s) in the strip < a < 1 are 
the points s = r for all square free positive integers k. 
Therefore, with the exception of those points, Eq. f) 13[) . 
together with Eq. (|15p . gives a representation for P(s) 
valid for Re(s) > 0. 

Let us use that C(^ s ) = z i where 2 is a complex vari- 
able. The function In z is analytic at every point of its 
Riemann surface. We have that dP(s)/ds = P'(s) can 
be computed for Re(s) > 0. From Eq. (fT3)l . one has 



k=l 



C(ks) 
k ((ks) 



(16) 



As mentioned above, the requisite for defining a regular- 
ized determinant is to compute the corresponding spec- 
tral zeta function and its derivative at s = 0, in order to 
obtain a well-defined functional integral associated with 
the scalar field. In the case of prime numbers spectrum, 
the functional W[0] is given by 

W[0] = iln m 2 P(.s)U=o + ^ P 00U=o- (17) 

Collecting our results we get the following pure non- 
existence theorem: 

The impossibility of extending the definition of the 
analytic function P(s) to the half-plane a < ; implies 
that free scalar quantum fields with the sequence of prime 
numbers spectrum cannot exist. 

It follows from the above theorem that, in the 
weak-coupling regime, the generating functional of the 
Schwinger functions cannot be defined either, since the 
generating functional for the free theory is meaning- 
less. As remarked before, such a result can be extended 
to higher spin fields, in particular spinor fields, with 
weak interactions. On the other hand, an intriguing is- 
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sue is whether this situation persists in non-perturbative 
regimes. 

Since the one-loop effective energy and the Casimir 
energy are closely related |25j . it is natural to discuss 
the vacuum energy associated with a theory in which 
prime numbers appear in the zero point-energy. Suppose 
a two-dimensional space-time and a box of size L. As- 
suming suitable boundary conditions, the vacuum energy 
is defined by (0|i?|0) oc J2p where the sum runs over all 
prime numbers. For the usual cases of regularization of 
ill-defined quantities, it is possible to prove the follow- 
ing statement: if the introduction of a exponential cut- 
off yields a analytical function with a pole in the origin, 
the analytic extension of the generalized zeta-function is 
equivalent to the application of a cut-off with the sub- 
traction of the singular part at the origin [26l - [28| . Once 
we accept the advantage of the zeta-function method over 
the cut-off method, we face a serious problem. Since the 
prime zeta function has an analytic continuation only in 
the strip < a < 1, one cannot obtain a finite renormal- 
ized vacuum energy for this particular situation. One can 
show that the same situation occurs in an interacting field 
theory. It is not possible to interpret the generating func- 
tional of Schwinger functions nor the one-loop effective 
action of the model in terms of well-defined mathematical 
objects. 

Similar conclusions can be drawn to quantum mechani- 
cal systems with n degrees of freedom since, as remarked, 



they can be derived from n-component scalar field the- 
ories in one dimension (time). Accordingly, the results 
obtained here put series restrictions on the formulation 
of quantum mechanical systems with a pure prime num- 
bers spectrum. On the other hand, suppose that there is 
a hypothetical quantum system realized in nature which 
has an energy spectrum given by the sequence of prime 
numbers. Hence from the above theorem one is forced to 
accept that there cannot be a free energy associated with 
such a system. 

It has been suggested in the literature a way to es- 
tablish the existence of a quantum system possessing the 
nontrivial zeros of the £(s) function as energy eigenval- 
ues by constructing one-dimensional potentials with a 
multifractal character. Whether or not this can be ex- 
tended to fields depends crucially on one's ability to set 
up such fractal potentials in the quantum field theory 
framework |29l - l3ll ]. This subject is under investigation 
by the authors. 
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